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Abstract 

Background: Myocardial T1 -mapping methods such as MOLLI use SSFP readout and are prone to 
frequency-dependent error in T1 -measurement. A significant error in T1 may result at relatively small off-resonance 
frequencies that are well within the region without banding artifacts. 

Methods: The sensitivity of T1 -estimates based on the SSFP based MOLLI sequence to errors in center frequency 
are calculated by means of a Bloch simulation and validated by phantom measurements. Typical off-resonance 
errors following local cardiac shimming are determined by field mapping at both 1.5 and 3.0T. In vivo examples 
demonstrate the artifactual appearance of T1 -maps in the presence of off-resonance variation. 

Results: Off-resonance varied 61.8 ± 15.5 Hz (mean ± SD, n = 18) across the heart at 1.5T and 125.0 ± 40.6 Hz 
(mean ± SD, n = 18) at 3.0T. For T1 = 1000 ms, the variation in T1 due to off-resonance variation was approximately 
20 ms at 62 Hz, and > 50 ms at 125 Hz. 

Conclusions: Regional variations due to the inability to completely shim the BO-field variation around the heart 
appear as regional variation in T1, which is artifactual. 
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Background 

The frequency response of steady state free precession 
(SSFP) sequences is well known and results in dark band 
artifacts [1-3]. Myocardial Tl-mapping methods such as 
the modified Look-Locker inversion recovery method 
(MOLLI) [4,5] use SSFP readout and are therefore prone 
to frequency dependent errors in Tl -measurements. The 
off-resonance behavior of SSFP and associated banding 
artifacts are often analyzed assuming the magnetization 
is at steady state as in continuous cine imaging. How- 
ever, the MOLLI imaging sequence uses single shot 
imaging with data acquired on the approach to steady 
state. As a result, the off-resonance response becomes 
dependent on the initial condition, which is in turn 
dependent on the inversion recovery time. This leads to 
a frequency dependent change in the apparent inversion 
recovery. It is not well appreciated that a significant 
error in Tl may result at relatively small off-resonance 
frequencies that are well within the region without 
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banding artifacts. As Tl -mapping and extra-cellular vol- 
ume (ECV) mapping based on Tl measurement are used 
to detect more subtle cardiomyopathies [6-14], small 
errors in Tl become more significant. Hence, even small 
variations in Tl introduced by off-resonance due to the 
inability to completely shim BO-field variation may lead 
to significant biases in measured Tl that may be falsely 
confused with real pathology. 

The readout excitation flip angle used in MOLLI is 
typically low (35°) compared to SSFP imaging as used 
for cine function imaging to reduce the influence of the 
readout on the inversion recovery. We hypothesize that 
reducing the flip angle improves Tl -measurement accur- 
acy with significant effect on the sensitivity to off- 
resonance. In this work, the sensitivity of Tl measurement 
to off-resonance frequency and readout excitation flip 
angle are quantified for specific MOLLI protocols. 



Methods 

Simulation 

A waveform level Bloch simulation of the MOLLI Tl- 
mapping method was implemented to study errors in 
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Tl -mapping, their dependencies and their sensitivities to 
various sequence and protocol design parameters. The 
simulation included the RF excitation pulse waveform 
and gradients in order to accurately model the variation 
in flip angle across the slice profile. Simulation results 
are provided for a specific MOLLI protocol which used 
a 5 (3s) 3 sampling scheme (similar to protocol used in 
[10,11]) which uses 2 inversions with images acquired 
for 5 heartbeats following the first inversion, followed by 
a 3 second recovery period, and images acquired for 3 
heartbeats following the second inversion. Tl fits are 
performed to the simulated magnetization using a 3- 
parameter model and a Look-Locker correction is applied 
to the apparent Tl to correct for the influence of readout 
[4,5,15]. The Tl-error due to the approximation used in 
the Look-Locker correction is the predominant error in 
this analysis [16]. Influence of magnetization transfer (MT) 
or imperfect inversion [17] are not considered here. 
Imaging parameters were: 256 x 144 matrix, FA = 35°, TR = 
2.8 ms, TImin = 105 ms, Tlshift = 80 ms. A partial Fourier 
acquisition in the phase encode dimension was used to re- 
duce the minimum achievable inversion time and reduced 
the overall shot duration. A partial Fourier factor of 7/8 
was used in the phase encode direction with 126 actual 
acquired lines (144 matrix size) plus 12 additional central 
lines for parallel imaging auto-calibration. The center of k- 
space was at 33 lines and there were 5 additional linear 
ramp flip angle pulses to reduce transient oscillations 
[18]. The Tl error was calculated for a range of Tl's 
from 200-1200 ms at a fixed T2 = 45 ms. The errors were 
also simulated for the specific values of phantom Tl and 
T2 used for comparison with measurements. 

Additional simulations were performed to assess the 
sensitivity of off-resonance errors to imaging protocol 
parameters. Sensitivity to off-resonance arises due to the 
influence of the SSFP readout during the approach to 
steady state. Both the total number of RF pulses (Ntotal) 
and the number of pulses to the center of k-space 
(Ncenter) will affect the Tl -mapping accuracy and sensi- 
tivity to off-resonance. The measured phantom and in- 
vivo protocols used Ncenter/Ntotal: 33/75, based on a 
matrix size with 144 phase encodes, Fourier factor = 7/8, 
and 12 extra central lines acquired for in-place auto- 
calibration of parallel imaging (factor 2). A range of proto- 
cols were tested and results are provided for a significantly 
shortened Ncenter/Ntotal: 16/48, based on a matrix size 
with 128 phase encodes, Fourier factor = 3/4, and no extra 
central lines acquired corresponding to a separate refer- 
ence line approach for auto-calibration. 

Phantom measurements 

Experimental validation was performed for a CuS04 
doped agar gel phantom using the specific MOLLI 
protocol at varying off-resonance frequencies in 10 Hz 



increments. Measured and simulated Tl are compared 
for a test tube phantom with Tl = 1197 ms and T2 = 47 
ms (1.5T Magnetom AERA, Siemens Medical Solutions, 
Erlangen, Germany). 

Field map measurements 

In-vivo measurements of off-resonance maps due to 
variation in BO-field were acquired to determine the typ- 
ical expected variation in frequency. Measurements were 
made in n = 18 subjects referred for CMR assessment of 
known or suspected heart disease at both 1.5 (age 46 ± 20) 
and 3.0T (age 36 ± 17) using a multi-echo GRE sequence. 
Field maps were estimated as a byproduct of water fat sep- 
arated image reconstruction [19]. A second order shim 
was performed over a local volume (box) encompassing 
the whole heart. Off-resonance frequency variation was 
measured in the left ventricle for a mid-ventricular, single 
short axis slice per subject. 

In-vivo imaging 

In-vivo Tl-maps were acquired for normal healthy vol- 
unteers to illustrate the apparent variation in measured 
Tl with off-resonance. Images were acquired at several 
center frequencies to demonstrate the sensitivity. Imaging 
was performed on 1.5T Magnetom AVANTO and AERA 
scanners and 3T Magnetom SKYRA scanner (Siemens 
Medical Solutions, Erlangen, Germany). This study was 
approved by the local Institutional Review Board of the 
National Heart, Lung, and Blood Institute and all subjects 
gave written informed consent to participate. 

Results 

Off-resonance variation 

Off-resonance frequency in the LV was measured at 
1.5 and 3.0T after shimming over a local heart vol- 
ume. At 1.5T, the mean off-resonance frequency in 
the LV myocardium was 20.3 ± 13.0 Hz. The maximum 
off-resonance in the LV was 61.8 ± 15.5 Hz (n = 18). 
At 3.0T, the mean off-resonance frequency in the LV 
myocardium was 15.4 ± 29.3Hz. The maximum off- 
resonance in the LV was 125.0 ± 40.6 Hz (n = 18). 

Simulation 

Figure 1 shows the response of SSFP readout vs fre- 
quency for various flip readout excitation angles shows 
the familiar dark bands spaced at 1/TR (T1/T2/TR = 1000/ 
45/2.8 ms). The left panel of Figure 1 shows the steady 
state response and the right panel shows the response for 
the transient approach to steady state corresponding to 
single shot imaging with initial full magnetization. The 
transient response has greater magnetization and different 
shape off-resonance response curves than the steady state 
response. The transient approach to steady state for single 
shot imaging during inversion recovery with different 
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off-resonance frequency (Hz) off-resonance frequency (Hz) 



Figure 1 Off-resonance response for a variety of flip angles for SSFP at steady state (left) and during the approach to steady state 
(right) for n = 33 RF pulses to k-space center with TR = 2.8 ms and 5 pulse linear run up. Note that the vertical axes have different scale 
since the steady state magnetization is reduced. 



initial magnetization for each inversion time image has an 
off-resonance response (Figure 2, left) which distorts the 
apparent Tl-recovery (Figure 2 right). 

The off-resonance error in Tl is shown in Figure 3 for 
a MOLLI 5(3)3 protocol with 35 degree flip angle for a 
range of tissue Tl values with T2 = 45 ms. Error in ms is 
displayed on the top panel of Figure 3 and in percent on 
the bottom panel. The Tl error and sensitivity to off- 
resonance frequency increase for higher values of Tl. At 
Tl=1000 ms, measured Tl varied by 10ms (1%) across ± 
50 Hz, and 20 ms (2.0%) across ± 75 Hz. This variation 
across frequency is in addition to the measurement error 
for on resonance tissue. 

Sensitivity to off- resonance for a nominal Tl=1000 ms 
was simulated for a modified protocol acquired with 
fewer RF pulses corresponding to smaller matrix size, 
higher partial Fourier factor, and eliminating in-place 
auto-calibration. The off-resonance sensitivity at 100 Hz, 
i.e., difference between on- and off-resonant Tl -estimates, 



was calculated to be 42 ms at 1000 ms (4.2%) using the 
protocol with Ncenter/Ntotal = 33/75, and reduced to 
27 ms (2.7%) using the protocol with Ncenter/Ntotal = 
16/48. The off-resonance sensitivity at 50 Hz varied from 
0.8% to 0.5%, for the cases 33/75 and 16/48 respectively. 

Phantom measurements 

Experimental data is in close agreement with the simula- 
tion (Figure 4) in the calculation of off-resonance de- 
pendence of Tl -estimation using MOLLI with SSFP 
readout. 

In-vivo examples 

An example study at 3T demonstrates the sensitivity to 
off-resonance. In this example, a local 2 nd order shim is 
used (Figure 5) to minimize off-resonance variation. 
Despite shimming, there is a residual uncompensated 
off-resonance variation at regions near the tissue-air 
interface that creates a detectable change in apparent Tl 



SSFP transient frequency response 



inversion recovery 




-100 -50 0 50 100 

Off-resonance frequency (Hz) 



1.5 2 2.5 3 

Inversion time (ms) 



Figure 2 SSFP off-resonance frequency response during the transient approach to steady state with inversion recovery (left) illustrates 
that the off-resonance response depends on the initial magnetization thereby influencing the apparent inversion recovery (right) plotted 
for 0 Hz (blue), 40 Hz (magenta), and 80 Hz (red). (MOLLI 5(3)3 protocol with FA = 35 deg, TR = 2.8, N = 33 to center, + 5 linear ramp). 
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Figure 3 Off-resonance error in Tl using MOLLI 5(3)3 protocol with 35 degree flip angle for a range of tissue Tl values with 12 = 45 ms, 
showing error in ms (top) and percent (bottom). The T1 error and sensitivity to off-resonance frequency increase for higher values of Tl. 
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Figure 4 Off-resonance error in Tl comparing phantom measurements (red) with Bloch simulation (blue) for a phantom with 
Tl = 1 197 ms and T2 = 47 ms using MOLLI 5(3s)3 protocol with a 35° flip angle. 
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Figure 5 Local 2 nd order shim prescription and field map measurement using multi-echo GRE approach at 3T which jointly estimates 
water, fat, and field map. Off-resonance variation is particularly pronounced at the tissue-air interface. 



(Figure 6) which is further accentuated with even small 
errors in center frequency. 

The example in Figure 7 illustrates the trade-off be- 
tween SNR and Tl -measurement bias by acquiring data 
at different flip angles. Simulation of Tl estimates vs. 
measurements using Tl = 1365 ms and T2 = 45 ms. 



Using a lower flip angle (FA) trades SNR (precision) for 
improved accuracy and reduced off-resonance sensitivity. 
Reducing the flip angle from 35° to 20° causes a reduction 
in SNR from 36 to 27 in the septum and 29 to 23 in the 
lateral wall. The Tl estimate in the septum is 1284 ms at 
35° and 1330 at 20° which agrees well with simulation 





Figure 6 Tl-maps acquired at different center frequencies using MOLLI at 3T. Despite the use of a 2 nd order shim in a local volume around 
the heart, off-resonant variation across the heart (Figure 5) leads to local variation in the apparent T1 as indicated by arrows. 
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Figure 7 T1 maps and raw inversion recovery images acquired at 3T using MOLLI 5(3s)3 protocol at various flip angles illustrating 
tradeoff between SNR and T1 -measurement bias. Simulation of T1 estimates vs measurements using T1 = 1365 ms and T2 = 45 ms. 



(Figure 7) using T1/T2 = 1365/45 ms. The Tl meas- 
urement error off-resonance is calculated by simulation 
(Figure 8) for this example. Tl bias error due to the effect 
of the readout on the inversion recovery curve is -48 ms 
and -81 ms for FA = 20° and 35°, respectively. Sensitivity 
to off-resonance over ± 100 Hz is 98 ms and 164 ms for 
FA = 20° and 35°, respectively. 



Discussion 

Tl underestimation 

It is generally assumed that shimming and center fre- 
quency adjustment are not significant problems at 1.5T. 
Maximum off-resonance averaged 61 ± 15.5 Hz across 
the LV (n = 18) leading to an apparent variation in Tl of 
0.5% for Tl = 400 ms, and > 1% for Tl = 1000 ms 



simulated off-resonance response vs flip angle 
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Figure 8 Using a lower flip angle (FA) trades SNR (precision) for improved accuracy and reduced off-resonance sensitivity. Simulation 
compares 20° and 35° readout flip angle for T1 = 1365 ms and T2 = 45 ms (representative values at 31). T1 bias error due to the effects of 
readout on the inversion curve is -48 Hz and -81 Hz for FA = 20° and 35°, respectively. Sensitivity to off-resonance over ± 100 Hz is 98 ms 
and 164 ms for FA = 20° and 35°, respectively. 
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(calculated using a Bloch simulation for the specific 
MOLLI protocol described). This represents the 
expected mean error in Tl (n = 18), however, off- 
resonance > 80 Hz was observed in 4 of 18 subjects, 
which resulted in more significant Tl errors (> 3%) 
which could be confused with regional variation in Tl 
due to pathology. In a large study of normal subjects 
[20], the Tl measurements were reported to have stand- 
ard deviation < 2% at 1.5T. However, that study mea- 
sured the mean per subject across the heart which tends 
to downplay regional variation due to local susceptibility. 
At 3.0T the off-resonance is typically greater, with mean 
off- resonance found to be -125 Hz across the LV, after 
local shimming. The variation in Tl at 3T from on- 
resonance with this protocol was calculated to be approxi- 
mately 3% at 125 Hz for Tl = 400 ms, and > 5% at 125 Hz 
for Tl = 1000 ms. Off- resonance >150 Hz was observed 
(7 of 18 subjects) leading to significant regional variation. 

In this study, the readout size was 256 sampled in the 
frequency direction in order to improve the spatial reso- 
lution to help mitigate partial volume effects at tissue 
interfaces, e.g, myocardium and blood. This readout 
length results in a TR of 2.8 ms. A shorter TR (and 
readout) would somewhat decrease the sensitivity to 
off-resonance. The off-resonance sensitivity is also 
dependent on the number of RF pulses to the center 
of k-space and total number of readouts, and can be 
reduced by decreasing the matrix size and greater 
partial Fourier. Modifying the protocol in this way 
may result in lower spatial resolution and degraded 
point spread function leading to Tl -errors caused by 
greater partial volume contamination of the myocar- 
dium by the adjacent blood pool. 

The use of separate reference line auto-calibration of 
parallel imaging eliminates 12 phase encodes per shot 
(6 to center of k-space) using the current widely used 
in-place calibration approach. For this reason, we have re- 
cently adopted a separate reference line approach which 
acquires the auto- calibration data at the end of the 
imaging acquisition such that the reference lines do not 
influence the initial magnetization. This has the additional 
benefit of reducing the overall shot time by 34 ms which 
is particularly important at higher heart rates to reduce 
cardiac motion blur. 

Reduction in excitation flip angle can reduce the off- 
resonance sensitivity at the cost of decreased SNR. Re- 
duction in flip angle also reduces the absolute bias error 
(i.e., on-resonance error). Using a FA = 35° the SNR is in 
the range 20-40 across the LV and the precision of Tl 
estimates is between 45 and 22 ms (SD measured on a 
per pixel basis), respectively, using the 5(3)3 protocol 
[21]. By reducing the flip angle from 35° to 20° the SNR 
is reduced by approximately 40% to the range of 15-30 
across the LV. With this SNR reduction, the precision of 



Tl estimates is 60 to 30 ms, respectively. The precision 
of the Tl estimated in a myocardial sector is improved 
several fold by averaging the voxels within the sector, 
however, per pixel precision is still important for 
detecting subtle regional variation. This trade-off might 
be worthwhile particularly at 3T where there is greater 
SNR to begin with and the magnitude of the off-resonance 
problem is expected to be greater. At some point, a 
FLASH based readout also becomes competitive but this 
requires further optimization, which is beyond the scope of 
this study. 

Off-resonance variation 

Off-resonance maps may be useful for improved con- 
fidence in measured Tl -values, although it would be 
difficult to correct the Tl values. Values for typical 
off-resonance variation were established in a relatively 
small patient population where the variation is pri- 
marily caused by tissue air interfaces. In subjects with 
sternal wires or devices the off-resonance variation 
may increase. In severe cases with rapid spatial variation 
of off-resonance, an additional issue may arise in cases 
where the subject has difficulty breath-holding. In this 
case, the possibility arises that the off-resonance variation 
over time will also affect the estimate of Tl, but this case 
has not been studied. 

Lastly, since it has been shown that the sensitivity to 
off-resonance is a function of flip angle, there is a com- 
plex interaction between the variation in off-resonance 
due to BO-inhomogeneity and the variation in actual 
flip angle due to Bl transmit inhomogeneity. The B0 
and B1+ field inhomogeneities are fairly independent, 
however sensitivity of Tl estimate to variation in flip 
angle has not been studied. 

Off-resonance sensitivity of ECV 

Errors in Tl due to off- resonance may also influence the 
calculation of ECV maps [10] derived from Tl-maps. A 
preliminary analysis has shown that the systematic bias 
errors in ECV due to off- resonance is quite low, on the 
order of 1% or less (in units of ECV percentage). Al- 
though the blood Tl is much longer, errors in ECV are 
dominated by the pre-contrast Tl. Interestingly, the Tl 
measurement blood is not sensitive to off-resonance as a 
result of flow. The source of error in the stationary myo- 
cardial tissue Tl estimate using the MOLLI method is 
due to the approximation of the so called Look-Locker 
correction (B/A-l) which arises due to modification of 
the apparent inversion recovery curve leading to a Tl* 
that is shorter than Tl. In the case of flowing blood 
where there are new spins at each beat, the inversion re- 
covery is not influenced by the previous readout and the 
Tl* from the 3-parameter exponential fit provides and 
an accurate estimate of Tl. 
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Conclusions 

The use of single shot SSFP imaging for Tl -mapping in 
the heart is accompanied by an off-resonance sensitivity 
of the Tl estimates. Significant Tl measurement errors 
may arise well within the SSFP "passband" in regions that 
do not experience dark band artifacts. The off-resonant 
behavior is due to the transient approach to steady state, 
which depends on the initial magnetization determined by 
the inversion recovery. The apparent Tl depends on off- 
resonance leading to a Tl-underestimation. Off resonance 
errors worsen with higher excitation flip angles. Regional 
variations due to the inability to completely shim the 
BO-field variation around the heart appear as regional 
variation in Tl, which is artifactual. 

Abbreviations 

CMR: Cardiovascular Magnetic Resonance; MOLLI: modified Look-Locker 
inversion recovery; Tl: Inversion time; ROI: Region-of-interest. 
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